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Abstract 

Quasi-one-dimensional stochastic Dirac operators with an odd number of channels, time 
reversal symmetry but otherwise efficiently coupled randomness are shown to have one 
conducting channel and absolutely continuous spectrum of multiplicity two. This follows by 
adapting the criteria of Guivarch-Raugi and Goldsheid-Margulis to the analysis of random 
products of matrices in the group S0*(2L), and then a version of Kotani theory for these 
operators. Absence of singular spectrum can be shown by adapting an argument of Jaksic- 
Last if the potential contains random Dirac peaks with absolutely continuous distribution. 



1 Introduction 

In this paper we consider a random family of Dirac operators H on the Hilbert space C^^) 
of square integrable functions with fibers of dimension L G N. It is of the form 

H = Jd + W + V,5.,, J = (^^ (1) 

where d is the space derivative, the potential W is a locally integrable function with values in 
the hermitian matrices Her(2L, C) of size 2L and Vj G Her(2L, C) are singular potentials at the 
points Xj G M (defined as usual by boundary conditions at Xj, see Section |2]). The potential 
W is a particular space-homogeneous random process described in detail below, and the Vj are 
independent and identically distributed. Both potentials are supposed to satisfy time reversal 
symmetry 

J*W{x) J = W{x) , J%J = Vj . (2) 

This means that jyV{x) and JVj are elements of the Lie algebra so*(2L) of the classical Lie group 
S0*(2L) given by those complex 2L x 2L matrices T satisfying T*JT = J and T*T = 1. Hence 
the Hamiltonian H is self-dual, namely J*HJ = H, and in the so-called symplectic symmetry 
class describing time-reversal invariant particles with odd spin. Apart from this symmetry, we 
suppose the coupling of the potential to be efficient. This is guaranteed if the distribution of 
J'Vj has an absolutely continuous component w.r.t. the volume measure on so*(2L), but can 
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also be satisfied by adequate choice of W if the V/s vanish. A more technical formulation of the 
(actually much weaker) couphng hypothesis is given below in Section [61 Our main new result is 
now: 

Theorem 1 Consider the random Dirac operator ([T]) with time reversal invariance ([2]) satisfying 
the Coupling Hypothesis on the randomness stated in Section |H 

(i) For even channel number L, the spectrum of H is almost surely singular. 

(ii) For odd channel number L, H has almost surely absolutely continuous spectrum of multiplicity 
2 on all ofR. If the distribution of the JVj is absolutely continuous on so*{2L) , the absolute- 
ly continuous spectrum of H is almost surely pure. 

Theorem [1] does not say anything about the singular spectrum in general {i.e. without the 
supplementary assumption on the distribution of the V/s), but we believe it to be always empty. 
It is crucial that L is odd, as discussed by several authors in the physics literature (please consult 
[EM] for a long list of relevant references). We believe that for even L the spectrum is almost 
surely pure-point, but did not try to prove this in detail (it should be possible by adapting the 
techniques of jKLSl IBouj ). The main difference between the odd and even case is that there are 
two vanishing Lyapunov exponents in the odd case and no vanishing Lyapunov exponent in the 
even case. This is related to Kramers' degeneracy and symplectic symmetry of the Lyapunov 
spectrum and is proved in Section [7l Based on this fact, the proof of Theorem [T] goes on by 
applying Kotani theory for Dirac operators as developed by Sun [Sunj along the lines of the 
work by Kotani and Simon |KSj . Even though most of the main identities in |Sunj are correct, 
it contains some errors which we felt necessary to correct here. Section [5] also generalizes the 
works [KSt ISunj to singular and complex-valued potentials. This extension of Kotani theory 
is non-trivial and crucial for two reasons: the Coupling Hypothesis cannot be satisfied for real 
potentials (see the arguements below) and the singular potentials are perturbations of finite rank. 
The latter leads to similar formulas for the Green functions as in rank one perturbation theory. 
Thus the last claim of the theorem can be proved by adapting the argument of Jaksic and Last 
|JL] (see Section [8]). Sections [2] to H] contain preparatory material some of which doesn't seem to 
have appeared in the literature and makes this work essentially self-contained. 

Let us put Theorem [1] in some perspective, both from a mathematical point of view and 
a physical one. Most quasi-one-dimensional discrete and continuous random Schrodinger op- 
erators exhibit Anderson localization, even though some peculiarities such as in the random 
polymer model may lead to non-trivial quantum diffusion [JSSj . The situation is different for 
first order differential operators. For example, consider h = 1 ^ td + v on L^(]R, C'^) where 
V G L°°(M, Her(L, C)) is an essentially bounded hermitian potential (which may be thought of 
as random). Then the initial value problem du = tvu, u{0) = 1, has a unique solution u = u{x), 
which lies in the unitary group U(L). Let us use it to define a unitary U on L^(M, C^) by 
{Uip){x) = u{x)%l){x). Then WhU = 1 ®id showing that h has absolutely continuous spectrum 
of multiplicity L for any potential v. In physical terms, the operator h can be thought off as 
an effective model for the chiral edge states of a quantum Hall system with edge conductivity 
L, and the above shows that the nature of the spectrum is conserved under perturbation by a 
potential, as is the Landauer conductivity which is equal to L (because U commutes with the 
position operator X on L^(M, C"^)). Note that the stability of the nature of the spectrum could 
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also be deduced from Mourre theory because i[h,X] = 1. For true edge states of a disordered 
magnetic operator on a half-plane, the proof of conservation of absolutely continuous spectrum 
|BPt IFGW] and the edge conductivity |KRSj is much more involved, but possible. 

Next let us explain why we believe that Mourre theory cannot be applied to the Dirac operator 
H because there is no natural conjugation operator. In fact, the only physically reasonable choice 
would be the spin current given by the time derivative of the self-adjoint observable A = iJX 
where is X is the position operator. However, i\A^ Jd + W] = 1 + XJ{W — J*WJ) is positive 
only if the time-reversal invariant potential W is real and thus JTW is in the Lie algebra so(2L). 
In this situation the Coupling Hypothesis is not satisfied and all Lyapunov exponents vanish. 
Theorem [T] is hence a much more delicate result than the one for h = id + v just described. 
We also find it to be a challenging problem to prove absolute continuity of the spectrum for 
half-plane models for which ([T]) is an effective description of the edge states. 

Next let us comment on the physical relevance of the Dirac Hamiltonian ([T]) with time reversal 
invariance ([2]). It is believed to be an effective model for so-called helical edge states in graphene 
sheets with a gap at the Dirac point (opened by spin-orbit coupling |EMj ). In such graphene 
sheets the number of edge channels with spin up and spin down is odd and hence these edge 
states are protected against localization. This is reflected by Theorem [H 
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in Cambridge. We also thank the Cambridge Philosophical Society for supporting the stay of 
Christian Sadel at the Newton Institute. This work was funded by the DFG. 

2 Weyl-Titchmarsh matrices 

This section introduces and analyzes Weyl-Titchmarsh matrices for a fixed non-random Dirac 
operator with point interactions. In part this is review (compare e.g. |HS] ) and therefore proofs 
are kept short, but results need to be written out if only to fix notations. Let § = (a;j)jgz be a 
discrete subset of M with no accumulation point and associate to each so-called singular point 
Xj a singular potential Vj G Her(2L, C). Furthermore let W be in the space Lj^^^(M, Her(2L, C)) 
of locally integrable functions with values in the hermitian matrices of size 2L. All this data 
encoded in = (W, (xj, Vj)j(zz), but in this and the next section u is fixed and hence suppressed 
in all notations. The time-reversal symmetry ([2]) is implemented only in Section [6l The first aim 
is to make mathematical sense out of H given in ([T]) as a self-adjoint operator on L^(M, C^'^). As 
usual, the singular potential is dealt with as a certain self-adjoint extension. Before going on, let 
us point out that most results of this paper also hold for the self-adjoint operator TZd + W where 
X t— > 7l{x) is bounded, invertible, and satisfies TZ* = —71 as well as dTZ = W* — W. In order to 
focus on the essential difficulties, we stick to the case TZ = J . 

Let W^^'^(M/S,C2^) be the Sobolev space of functions L2(M/§,C2^) with square-integrable 
first distributional derivative. Note that these functions ip are continuous away from S and have 
left and right limit values il){x±) = lim^ioip{x ±e) for all x G M. First we consider the restriction 
Hq = H\'D{Ho) to the domain 

V{Hq) = {V e iy^'^(M/S,C2^) I ^(x+) = %Ij{x-) = for x G § } . 
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Then the domain of the adjoint is T>{Hq) = H^^'^(R/§, C^'^). The proof of the following result is 
adapted from |LMj . 

Proposition 1 Forip,(l) E 'D(Hq), one has 

where the scalar product on the l.h.s. is in L^(R, C^-^) and those on the r.h.s. in C^^. 

Proof. Let Xn e C°°(M, [0, 1]) with Xn I [-«,«] = 1, Xn\[-2n,2nY = and x'n = dXn < f for some 
constant C. For any G V{H^) set <pn = Xn(t)- Then 0„ ^ and ifo*0n H^cj) in ^^(R, C^-^). 
Therefore one can calculate as follows: 



where we used the local integrability of W. This directly implies the proposition. □ 

If § is empty, then the r.h.s. of ([3]) vanishes and this shows that Hq is self-adjoint with 
domain iy^'^(R, C^'^). In the terminology of Weyl theory described below, this means that H is 
in the limit point case for any locally integrable potential W. This fact also follows from Weyl 
theory (more precisely, the bound fll2p below) without reference to Proposition [H If § is not 
empty, then Hq has non-trivial deficiency spaces (which are infinite dimensional if and only if 
S is infinite). Beneath all the self-adjoint extensions of Hq we are interested in those given by 
local boundary conditions, namely those not mixing the deficiency spaces corresponding to each 
of the terms on the r.h.s. of ([3]). Within the class of local boundary conditions we will choose 
the ones obtained by formally approximating the singular potential Vj5xj (this will be explained 
below), namely we consider the domain 

V{H) = G Vr^'2(R/§,C2^) I ^l){xj+) = e^^^-^ixj-) for j G N } . (4) 

Then H = Hq\x>{h) clearly is an extension of Hq and the identity (e^^^)*i7e"^^J = J' replaced in 
([2]) shows that it is self-adjoint. 

Now that the operator H is well-defined, let us introduce the transfer matrices (or fundamental 
solutions) T^{x,y) GMat(2L x 2L,C), x > ?/ G R, at a complex energy 2; G C as the unique 
solutions of 

{H-z)ri.,y) = 0, T^{y,y) = Ul , (5) 

which are right-continuous in x and in y (for x > y) and for which x ^ T^{x,y) is in T>{H). 
(Recall that a function is left-continuous if f{x—) = f{x) for all x and right-continuous if f{x+) = 
f{x) for all X.) For x < y, we set T^{x,y) = T^{y,x)^^. At xj G § the transfer matrices then 
satisfy T^{xj,Xj—) = e"^^^. The general composition rule reads for x,u,y E M 

r^{x,y) = r{x,u)r{u,y) . (6) 



= lim {HQipn I 0m) - (V'n | -f^o0m) 
n,m—> 00 
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For later convenience we also set T^{x) = T^{x,0). Now let us briefly sketch in which sense 
the boundary conditions in (jlj) are natural. Indeed, if Xn € C^(M, M) converges weakly to Sxj 
and T^{x,x') is the transfer matrix with potential VjXn, then taking the limit n — > oo first, one 
formally verifies T^{xj, Xj—) = e'^^^ which is precisely the jump condition above. Next comes 
the basic but crucial Wronskian identity for the transfer matrices. 

Lemma 1 For a < b and 2;, G C, 

r'{h-YjT^{h-) - r\aYjr^{a) = (c-^) / dxr{xyT^{x). (7) 

J a 

Proof. Denote the points in § U (a, 6) by xi,...,xn and set xq = a and xn+i = b. Then 
X I— > T^{x) is differentiable away from these points. Thus, using the local integrability of W, 



{(-z) I dxr\xyr^{x) = ^ /" '^'dx [r"(x)*((c- w)r^(x)) - ((2- w)r"(a;))*r^(x)] 

N 

= E [r{x,^i-rjr'^{x,+i-) - r\x,+yjr^{x,+) 



j=0 "^J 
N 

j=0 

where the second equality follows from the differential equation ([5]) and the fundamental theo- 
rem. Replacing T^(xj+) = e'^^^T'^{xj—) and using (e'^^j)*j7'e^^^ = J', one sees that only the 
boundary terms remain and thus the lemma follows. □ 

Next let us consider the restrictions of H to M+ = (0, oo) and ]R_ = {—oo, 0) given by 
H± = if|2,2(]R^ c2L). These operators are not self-adjoint because the same calculation as above 
shows 

{H*^^P\<P) - {^\H*^(p) = ±i^iO±rj<P{x±) , (8) 

for e V{Hl) = e Vri'2(R±/§,C2^) I tl){xj+) = e^^^ip{xj-) for j G N }. This shows 
that the self-adjoint boundary conditions for H± are precisely given by the set of hermitian 
Lagrangian planes, namely = {$ G Mat(2L x L, C) | rank($) = L,^*J^ = 0}/ ~ where 
$ ~ $' $ = $'c for c G GL(L,C). For one such plane $ G hi, the associated self-adjoint 
operator will be denoted by H± is,. It is well-known (see e.g. [SBl] for a short proof) that hi is 
diffeomorphic to the unitary group U(L). Thus the deficiency spaces A^^ = keT{H^ — z) of H± 
are L-dimensional. 

For any analytic function g we denote its complex derivative by d^g = g. 

Theorem 2 For '^m{z) ^ there exist unique so-called Weyl-Titchmarsh matrices G 
Mat(L X L,C) such that ker(if^ — z) is spanned by the column vectors of 

= rix) ( ^]^, ^ . (9) 



{Here the column vectors of are considered as elements of L'^{W±,C'^^), but below is 
also used for all x G M.) They are analytic in C/M and satisfy the Herglotz property 

= [ dx ^u^r^u^) > , (10) 
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where '^m{Z) = — Z) for any operator Z , as well as 

{Miy = Ml, Ml = [ dx ^lixy^iix) . 

Proof. Let us consider the case of the sign + and '^m{z) > 0. It was argued above that the 
dimension of ker(if^ — z) is L. As every solution of H+ip = zip is of the form ijj[x) = T^{x)v 
for some vector v G C^^, it follows that there are L x L matrices a and (5 such that the column 
vectors of 

;w ) - ^'<^) ( ; 

span ker(ff^ ~ z)- As these vectors are, in particular, square integrable, replacing them twice in 
the Wronski identity ([7]) with b = oo and a = shows that 

POO 

i{l3*a-Q*(3) = 2^m{z) / (a(a;)*a(x) + > 0. 

Jo 

From this it follows that both a and (3 are invertible because for a vector v in the kernel of a or 
(3 one would have v*{(3*a — a*f3)v = 0. Therefore one can set M^ = f3a~^ and this also leads to 
the formula (ITU]) . The identity (M|.)* = M|. follows by replacing C = ^ and a = 0, 6 = cxd in the 
Wronski identity (JTj). Finally, let us check the analyticity of and derive the formula for its 
derivative. Again the Wronski identity with a = and b = oo shows for z ^ ( that 

^ A*^C/ A Mi -Ml 



Note that the integrand on the l.h.s. is square integrable also in the limit ( ^ z {at least for 
z G C/M), so that is indeed holomorphic and the formula for the derivative follows. The 
proofs for Ml are similar. Let us point out though that due to our definitions the jump 6"^^° is 
relevant for Ml if xq = G S. This is of some importance below. □ 

As a short aside, let us sketch how the modeling of the singular potential in ([1]) by the jump 
conditions in (jlj) fits with the theory of extensions by von Neumann. For this purpose, let us 
add the singular potential V = Vq at Xq = to the operator H. Let Hq be the restriction of H to 
V{Hq) = {'?/' G Vi^H) I ^/'(0+) = TpiS^—) = 0}. Due to Theorem [2] the deficiency spaces are both 
2L-dimensional and given by ker(i7o - z) = ^^C^ © ^^C^, where 

1 

n(^) = x{±x > 0) r{x) ( ) {Ml - {Mir: 

and X is the indicator function. These are partial isometrics : C'^ N^, namely \E'^(\1'^)* 
is the projection on and (\E'^)*\E'^ = 1^. Now the unitaries from ker(i7o — z) to ker(ifo — 'z) 
parameterize the self-adjoint extensions of Hq. Using the partial isometrics, these unitaries are 
precisely given by \I/1)?7(\E'+, ^E'l)* where U runs through the unitary group U(2L). It is 
now a matter of calculation to check that 

U = [(vl/^(0+),0)-e^^(0,vI/!(0-))]-' [(xl/;(0+),0)-e^^(0,M/i(0-))] , (11) 
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is well-defined {i.e. the inverse exists), is unitary and gives exactly the self-adjoint extension 
given by the jump condition ijj{0+) = e'^^'?/'(0— ). Hence every local boundary condition in (jl]) is 
an extension within the local 2L-dimensional deficiency spaces in the sense of von Neumann. On 
the other hand, there are local von Neumann extensions which are not given by jump conditions 
(for example, those which do not couple left and right). 

Even though it was already shown above that H is always self-adjoint (so that one is always 
in the limit point case), we now describe the Weyl theory because it gives quantitative estimates 
for the Weyl- Tit chmarsh matrices needed below. We closely stick to the notations of our prior 
work |SB2] along the lines of which also the proofs of the results below can be given (even though 
there are definitely older references such as |HSj for some of them). The basic idea is to study the 
restriction of the operator to L'^{{0, x), C'^^) and to analyze which initial conditions at lead 
to solutions satisfying any self-adjoint boundary conditions at x (there is an analogous treatment 
for H_). If an adequate chart for these initial conditions is used they have the geometric structure 
of a matrix circle in the upper half-plane, called the Weyl surface. As x increases, this circle 
shrinks in a nested manner. In the so-called limit point case that one always encounters for 
the Dirac operators, it shrinks to a single point in the limit x oo identified with the initial 
condition of ([9]) specified by the Weyl-Titchmarsh matrix M^. This fact reflects that there is 
no need to fix a boundary conditions at infinity in this case (the L^-condition takes care of it) 
because H is already self-adjoint. 

Now comes the more technical description. Let Gl be the Grassmannian of L-dimensional 
planes in C^^. The chart on used is the stereographic projection vr sending an 2L x L matrix 
(^) representing the plane to aP~^ G Mat(L, C). It is defined on full measure subset G'l'^ C 
on which the inverse of (3 exists. Then the Weyl surface at x 7^ is defined by 



dW{x) = -7r({$ G Gl |T"(x)<I) G Ll }) = \-M 



T\x) ( ^ ) G L, 

where the equality follows by showing that every plane $ in the first set is of the form in the 
second one jSB2t Prop. 7]. Now it is useful to rewrite the condition T^(x) $ G in terms of 
the quadratic form 

Q\x) = - r\xYjr\x) , 

namely 9W^(a;) = — vr ({$ G Gl | $ isotropic for Q^(a;) }). The definition of Q^(x) shows that 
Q^(a;+) = Q'ix—) also for a; G S so that Q^{x) is continuous and thus 9W^(x-|-) = c}W^(x— ). 
Item (i) and (ii) of the following properties of Q^{x) follow from the definition and the Wronskian 
identity, while (iii) can be checked as in |SB2j once one has verified that T^x)'^ = J*T%x)J. 

Proposition 2 The quadratic form Q^{x) satisfies: 

(i) Q'{x) = '-J + 2 Qm{z) (T-( . )\T% . ))LHio,.),c^^) 

(ii) ^m{z)dQ'{x) > 

(iii) Q'{x)-^ = J*Q~ix)J 

Now the radial and center operator are defined by 

r/ 1 \* / 1 

R'ix 



S'{x) = R\x) [ J ) Q\x) ( I 
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Both R^{x) and S^{x) are continuous in x (apart from the singularity at x = 0). It follows from 
item (i) of Proposition [2] that R^{x) > and —R^{x) > for Qm{z) > 0, and item (ii) implies 
dR^{x) < 0. The terms radial and center operator are justified by the following result which can 
be checked by the same calculation as in [SB2j . It is the basic fact of Weyl theory. Let the matrix 
upper half-plane be defined as the set of matrices Z G Mat(L, C) satisfying 53m(Z) > 0. 



Theorem 3 Let '^m{z) > 0. Then 

dW'{x) = I + i?"(x)5?7(-i?^(x))5 



U*U 



C V, 



If now the open and closed Weyl disc W^(x) and W^(x) are defined by this formula with U 
running through the set defined by U*U < 1 and U*U < 1 instead of the unitary group U(L), 
then the Weyl surfaces are strictly nested in the sense that for x > x' > or x < x' < 



[x) C W"(x') , dW{x') n W^(x) 



This theorem can also be used to prove the uniqueness of instead of the above argument 
based on ([H]), that is, basically the calculation in the proof of Proposition [H Indeed, along the 
lines of Proposition 11 of [SB2j one can prove that there exists a constant c such that 

11^^(^)11 < I I J . (12) 
\x\ ^m(zj^ 

This implies that H± is in the limit point case in the literal sense and that one furthermore 
has — (M^)~^ = lim^-^-i-oo 5*^(0;). We need the following consequence for our purposes below. It 
replaces perturbative arguments in [KSl ISun] and hence the bounds below hold under the more 
natural assumptions that W is locally integrable. For Schrodinger operators a similar reasoning 
applies if they are supposed to be in the limit point case. 



Corollary 1 There are constants Ci,C2 depending only on z and the Lj^^-norm ofW such that 

< Ci , - < ^ \ f < 02 . 

02 ^m[z) 

Proof. At x = the radial operator is infinite in the sense that R^{0)~^ = 0. As 

d{R^{x)-') = Qm{z) I ) rixyrix) ( J ) 

is equal to '^m{z)l > for x = and is continuous in x (even differentiable), it follows that 
R^{x)~^ > for some x > 0. Hence ||i?^(x)|| < 00 and the Weyl disc W^(x) is compact 
and strictly contained in the upper half-plane U^,. Furthermore by Theorem [3] the limit point 
— (M|.)~^ is an element of W^(x). As Z 1— >• —Z^^ maps compact sets of Ul to compact sets of 
Ul the proof is complete. □ 
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3 Green's function and spectral analysis 



This section deals with the Green function and spectral theory of the self-adjoint operator ([T]) 
defined by (jlj). We always assume that xq = G S, set V = Vo and denote the operator with 
singular potential V by ify (hoping that the reader can distinct Hq with V = from the Hq in 
the last section). 

Proposition 3 Let Qm{z) ^ and M'^, T^{x) and he associated to Hq {this only leads to 
changes for x < and the sign —). The resolvent {Hq — z)^^ is an integral operator with kernel 

G^Q{x,y) = <^l{x){-Ml-Ml)-'<^l{yr , (13) 

where the upper and lower signs are taken if x < y and x > y respectively. Furthermore, for a 
Lagrangian plane $ = (1 7)*, the resolvent (i^+,<i> — z)^^ is an integral operator with kernel 

( T^(x) <I>(-M; + 7)-i $;(?/)* , x<y, 

g;_^(x,i/) = <^ 

[ <^l{x) (-M; + 7)-i $* T^{y)* , x> y . 

Proof. Let Gq be defined by the formula in the theorem. Using (M^)* = one readily verifies 
that for all x G M, 

lim\G'Q{x + e,x)-G'Q{x-e,x)] = r{x)jr{xy = J, (14) 

ej,0 

where the last equality follows by taking the inverse of T^{x)*J'T^{x) = J', which is the Wron- 
skian identity ([7]) with ( = z, a = and b = x. Therefore setting iIj{x) = J dy GQ{x,y)(j){y) 
for a smooth function (p G L^(]R, C^^), the definition ([5]) of the transfer matrices implies that 
{Hq — z)iIj = (j) because dsgn= 2Sq if sgn is the sign function and 6x is a Dirac peak at x. Hence 
Gq is indeed the desired integral kernel. The formula for the half-sided operator is verified in a 
similar manner. □ 

From Proposition [3l ffTTl) and the general Krein formula for resolvents of self-adjoint extensions 
one could now deduce an explicit formula for the integral kernel G\;{x,y) of Hy. Then lengthy 
algebraic calculations lead to Proposition H] below, but we can also deduce it more directly based 
on the following idea. Both functions x ^ Gy{x,y) and y 1— Gy{x,y)* = Gy{y,x) are in the 
domain V{Hi;) and satisfy respectively (ify — z)G^{., y) = 6y and {H\; — z)G^{x, .) = Sx- Away 
from xq = 0, the domain of V{Hq) and the identities for Hq are the same. Thus a good Ansatz is 

GUx,y) = Gl{x,y) + Gg(x,0+)/CGg(0-,y) , 

with a matrix /C to be determined. The jump condition G\,{Q+,y) = e'^^G^{0—,y) gives for 

Gg(0,y) + Gg(0+,0)/CGg(0,y) = e^^[G^(0,y) + Gg(0-, 0) /CGg(0, y)] . 

Now let us take the difference of this equation for y = 0+ and y = 0—. Because Gq(0-|-,0) — 
G'{0-, 0) = J hj one obtains 

J + Gg(0+,0)/CJ = e^^[J + G'q{0-,0)}CJ] 

This equation can formally be solved for /C, leading to the following formula. 
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Proposition 4 Let '^m{z) ^ 0. The resolvent {Hy — z) ^ is an integral operator with kernel 
GUx,y) = G^,{x,y) + Gg(x,0) [e^^Gg(0-,0) - Gg(0+,0)]"' (1 - e^^) Gg(0, y) . (15) 



Proof. It remains to check that the appearing inverse is indeed well-defined. Due to (fT3i) . 
there exist two L-dimensional planes $± with ±7r($-|-) G Vl such that Go(0— ,0) = $+$1 and 
Go(0-|-, 0) = Now we claim that for any hermitian symplectic T satisfying T*J'T = J', 

thus in particular T = e^^^, one has T$_|_C^ fl $_C'^ = {0}. This implies as desired that 
T$+$l — is invertible. To prove the claim we first note that 7r(T$_|_) G (as the Mobius 

transformation with a hermitian symplectic matrix sends Vl to Vl) so that it is sufficient to 
consider the case T = 1. Now let = for some v,w & C^. Set a± = (1 0)$-|- and 

P± = (0 1)$±) both of which are known to be invertible. Then a+v = a^w and = 
Thus v = P^^p -W so that aP^'^P-W = a^w. Therefore u = P-W satisfies a+Pj^^u = a^P_^u and 
thus «*??($+)« = «*??($_ )m. By hypothesis this implies u = and consequently w = v = 0. □ 

Before going on let us discuss the discontinuities of Gy{x,y) in the vicinity of the point 
{x,y) = (0,0) (any other singular point can be analyzed similarly). Because x G^{x,y) and 
y Gy(x, y)* = G^iy^x) are in the domain V{H\;), the singular potential leads to jumps on 
the lines x = and y = 0. According to (IT^ there is furthermore a jump by on the diagonal 
X = y. Away from these 3 lines crossing at the origin, Gy{x,y) is continuous. Hence there are 
6 directional limits as {x,y) (0,0). Enumerate them by Gi, . . . , Gg in a clockwise direction 
starting with Gi = \im^iQGl;{e, 2e). Setting T = e'^^ one then has 

G2 = G\ + JT" , = G2{T ^)* , (j4 = T , G5 = G4 — JT" , Gg = G^T* , Gi = TGq . 



Note that these relations are indeed cyclic because T*J'T + J . By (1131) each of the Gj has rank 
L. The following proposition shows that, however, an adequate linear combination is a Herglotz 
function and, in particular, of full rank 2L. 



Proposition 5 Let us define the averaged Green matrix 
G\,{x) = lim 



40 



^ G'y{x + e, X - e) + ^ G'y{x -e,x + e) + ^ G'y{x + e,x + 2e) 



-G^{x + 2e,x + e) + - G^{x -e,x -2e) + - G^{x - 2e, x - e) 



Then z eUi ^ G^{x) = (Gy(x))* G Mat(2L, C) is a Herglotz function for any x G R/S and has 
non-negative imaginary part for x E E>. It satisfies 

Qm{GUO)) = J (1 + e^"") '^miGUO+)) (1 + e^)* . (16) 

Proof. Let us note that for x ^ S the definition of the averaged Green matrix reduces to 
Gy{x) = |(Gy(x-|-,x) + Gy(x— ,x)). For sake of notational simplicity, let us focus on the case 

X = with V 7^ modeling x G S. With the above notations, then by definition Gy(0) = 
|(Gi + G2 + 2G3 + G4 + G5 + 2Gq) which is a weighing of the Gj according to the area of the 



10 



corresponding octant or quadrant. Now let z = E + te with e > and consider the positive 
operator Sm((i/v - z)'^) = e((ifv - + e^)-^ For any Lp G ^^(IR, C^-^), one thus has 

< {^\'^m{{Hy - zr')\^) = Jdx J dy ^{xr{GUx,y) -GUy,xr)^{y) . 

Now let Xk £ C'^(M) be a positive approximate unit, that is w-hnifc^oo Xfc = ^o- For any 
function / : = C ^ C having the directional limits f{9) = limrj^o /('^e*^), it follows that 
Jdx J dyxk{x)xk{y) f{x,y) converges to / = /g'" |f f{9). Hence, for ipk = XkV with v G C^^, 

< lim{^k\^m{{Hv - zy')\^k) = v* - Gm*) v ■ 

This proves that the imaginary part is non-negative. The Herglotz property for ^ S, namely 
that the imaginary part is positive, follows from the concrete formula 

^ / (-M+-M_)-i (-M+-M_)-i(M+-M_) \ 

V (M+-M_)(-M+-M_)-i (M-i + M-i)-i J ^ ^ 

following from Proposition [31 and the Herglotz property of M± by the Liouville theorem. As 
the singular points are discrete, there is an interval (0, e) not containing any. Hence Gy(0+) = 
|(G'i + G2)- It is now a matter of an algebraic calculation to verify the second formula. □ 

As for any Herglotz function with sufficient decay properties such as G^{x), there is associated 
a matrix valued measure fix on M and a self-adjoint matrix = A* independent of z (see |GT] 
for a review and properties) such that 



Because 

GUx) = T^ix,y)GUy)rix,yr 

for x,y ^ §> and T^{x,y) is analytic and invertible, the measures fix, x ^ S, all define the 
same measure class. According to ffTBl) . the measure /iq associated to Gy(0) is also in the same 
measure class as long as —1 is not in the spectrum of e^^^. We skip the proof of the following 
result, showing in which sense fix can rightfully be called a spectral measure of Hy (see |KSj ). 

Proposition 6 Let ip,4> E L^(M, C^"^) and f G Co(]R). Then, whenever fix is in the almost sure 
measure class, 

{^\f{Hv)\<p) = J^f{E) (^j dyr'^{y,xr^{y)^ fix{dE) (^j dyT'^{y,xr<P{y) 
and the functions of E in the parenthesis are in L'^{R,fix). 
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The arguments in Section [8] will be based on the following perturbative formula for the 
averaged Green matrix w.r.t. the finite rank perturbation given by the singular potential VSq. 
For notational convenience let us set = G^{0). Furthermore let us introduce the Cayley 
transform of V by 

V = 2:r(e^^ + l)-i(e^^-l) , (18) 

whenever the inverse is well-defined. One readily checks that V* = V and that J*V^J = V if 
J-YJ = V. 

Proposition 7 The averaged Green matrix satisfies (even ifV is not well-defined) 

^1 



o) 



and 



-1 



V 



1 + VG' 



(19) 
(20) 



Proof. Let us apply the averaging procedure of Proposition [5] to (ITS!) . This gives 

PiZ I PiZ Pz Pz /-I I 



where /C = [e^^Gg(0-, 0) - G^(0+, 0)] (1 - e^^) as before. Because both G^ and G§ are 
invertible, it follows that also (l + /C Gq) is invertible. Hence 



Gy = Gq + Gy ( 1 + /C Gq) ^ /C Gq = (Gq) ^ — (l + /C Gq) ^ K, 



(21) 



Using Gq(0±, 0) = Gq ± I J", one readily checks (l + /C Gg) ^/C = —V completing the proof of 
(fT9|). That of ([2Q]) is straightforward. □ 



4 Stochastic Dirac operators 

In this section we introduce stochastic Dirac operators and state a few of their elementary prop- 
erties, then introduce the random Dirac operators and give a precise statement of the main 
coupling hypothesis needed in Theorem [H Let be given a compact dynamical system (fi,P,T) 
where T is a continuous M-action on the compact space f2 w.r.t. which the probability measure 
P is supposed to be ergodic. Then {H{ijj))^^q is called a family of stochastic Dirac operators 
if each H{uj) is of the form ([1]) and the map uj & VL ^ H[uj) is strongly continuous in the 
resolvent sense and covariant, that is, if Ux denotes the right shift by x on L^(M, C^'^), then 
Ux{H{ijj) — z)^^Ul = {H{TxUj) — z)~^ . Each point G is thought of as a configuration, in- 
corporating the positions S and values {yx)x& of the singular potential as well as the potential 
W. Thus S is an R-ergodic point process. Its density is denoted by p§. The locally integrable 
potential associated to a given configuration oj is then >Va,(x) = >V(T_a,ci;), a; G M, where the 
W is a matrix-valued function on VL. Hence we suppose this function W to be locally integrable 
along orbits with a uniform bound on the L^-norm over unit intervals. 
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Now all objects such as transfer matrices, Weyl-Titchmarsh matrices and Green matrices 
analyzed in the sections above depend on cu; however, in the notations this will be made explicit 
by a supplementary argument only if necessary. Let us introduce some notations for the L x L 
matrix entries of the potential: 

\R* Q J ' \C D 

All these objects are random and for V = V^, x e §, the entries are also denoted A^, Bx,Cx, D^. 
As the matrix is in SP(2L, C), it is well-known that the inverse in the definition of the 
Mobius transformation 



(^^ d)-^ - {AZ + B){CZ + D)-\ 



exists whenever Z is in the upper or lower half-plane, i.e. ±'^m{Z) > 0. If then W = ■ Z = 
[AZ + B){CZ + D)~^, also W is in the upper or lower half-plane respectively and one has 
Z = (e^^)-^ -W = {D*- B*W){-C* + A*Z)-\ Now we can collect a few first properties of the 
transfer matrices and the Weyl-Titchmarsh matrices. 

Lemma 2 Let '^m{z) ^ 0, set 

(Sm) = ^'(^'^)(±4h) = *i'"'"'- '''' 

(i) The transfer matrices satisfy the cocycle equation 

r{x + y,u) = r{x,T_yu)r{y,u) , T%Q,u) = 1 . 

(ii) One has 



Pl{x + y,uj) ) \ Pl{x,T_yUj) 
In particular, a^{x,u!) is a cocycle: 



al{y,u;) . 



a^^{x + y,uj) = a^{x,T_yUj)a^^{y,uj) , a±{0,uj) = 1. 
{in) Ml{T^,uj) = ±pl{x,ij)ai{x,u)-^ . 

(iv) The map x i— M^{Txu) is differentiahle away from S. It is left- continuous and for —x G §, 

(v) The maps y h- >• Q;j_(a;, T^^a;) and y i— >• (3j.{x,TyU!) are left- continuous. For —y e §, 
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(vi) The map x G M+ ^ q;j_(x, tu) is right- continuous. If x G S, 

al{x,u) = {A^±B^Ml{T^.,+uj))al{x-,u) = {DIt BlM^T.^u))-' a'^{x-,uj) . 

(vii) d,ai{x,u) = [-R{T_,uy T {Q{T_,u) - ^)M|(T_,cu)] cu) /or x ^ § 

(viii) The following Ricatti equation holds for x ^ E> 

Proof, (i), (ii) and (iii) follow immediately from ([6]) and ( l22l) . It is clearly sufficient to analyze 
the directional continuity in (iv) and (v) for the case x = G S. Let e > 0. Using the composition 
rule for transfer matrices and their translation property 

T^x + e,y + e,T,uj) = T'{x,y,uj) , 

one deduces 

r{x,cu) = r{x + e,x,uj)-^r{x,T^,uj)r{e,0,uj). 

Taking the limit e | gives T^{x, u) = T^{x, Tq_uj) which implies M^{Tq_uj) = M^{uj). Similarly, 
the limit e | of 

T'{x,uj) = T%x,x-e,uj)T'{x,T,uj)r{0,-e,uj)~\ 

leads to 

T'{x,u) = e^^-rix,To+co){e^^'')-' . 
As the jump at x does not effect the square-integrability in ([9]), this implies that 

( ±Ml{u) ) ^ = ( ±M|(To+^) 

for some invertible L x L matrix A^. The upper entry implies that N = {D^ — Bq{±M^{uj)))^^ , 
the lower one 

± MliTo+uj) = i-C; ± AlM^iuj)) (D; t B^M^iuj))-' . (23) 

This is precisely the equation claimed in (iv) in the case x = 0. (v) follows from fl22l) and the 
last 4 identities. For (vi) we use T^{x,u!) = e^^'^T^ix—.uj) for a; > 0, giving 

^l{x,u) = e^^^ ^l{x-,u) = e^^^ ±M|(r_(,_)u;) ) «±(^-'^) ' 

where (iii) was used in the second equality. The upper entry of this identity gives the ffist equality 
of (vi). The second one follows by replacing ([23]) and using A^Dl-B^Cl = 1 and A^B* = -Bx^^- 
The following calculation gives (vii): 

d,al{x,u) = (10)5.r^(x,a;) ^^^^^^ ^ 

= (0 1)(^->V(T_..)) (±M|(ko.)) 
14 



Finally, 



so taking (vii) into account gives (viii). □ 

The fact that (l23l) is a Mobius transformation with a matrix out of SP(2L, C) has a number 
of consequences which we regroup for later use. 

Corollary 2 Let x e S and set M| = M|(T_,cj), M|(+) = M|(T_,+cj) and V = V,. T/ien 

i) ±M|(+) = {-C* ± A* M^) {D* T B* Ml)-^ 

ii) ±M|(+) = (±M|5 - D)-i (C ^ M|A) 

iii) M;(+) + (+) = (£> - M^B)-^ (M^ + M^) (D* + B*A'F)-^ 

= {D + MiBy^ + Mi) {D* - B*Ml)-^ 

iv) ±M| = (A±5M|(+)) (C±DM|(+))-i 

v) A ± 5M|(+) = (D* T 5*M^)-i 

vi) ^>m(M|(+)) = (DTM|5)-i$5m(M|) ((DtM|5)-i)* 

vii) S5m(M|(+)) = ((D* ^ fi*Af|)-i)* $5m(M|) (£)* ^ 5*M|)-i 

viii) M^(+)-Mf(+) = {MlB~Dy^M^{B*M^-D*)-^-{M^B + D)-^M^{D* + B*M^)-^ 

Proof. All this follows by short calculations using e.g. the Appendix of |SB2] and the identities 
AB* = B A*, CD* = DC* and AD* - BC* = 1. □ 

Now let us recall the definition of the Lyapunov exponents and state some of their properties. 
Because T^(x, u) is a cocycle by Lemma[2], Osceledec's theorem (see [KSj for a concise statement) 
associates 2L Lyapunov exponents at +cxo and — oo which will respectively be denoted by 7f > 
••• > 12L ^^"^ 7i — ••• — 72L- Similarly, a^{x,uj) are other cocycles of L x L matrices, so 
again each has L Lyapunov exponents at +oo and —oo denoted by 7^'^ > • • • > 7x,'^ and 
7^^'^ > . . . > ■ Part of the following proposition is copied from [KB] (even though the 
definition of 7^'"'' differs by a sign). 

Proposition 8 The various Lyapunov exponents satisfy: 

i) 7f = -I2L-1+1 for 1 = 1,..., 2L 

ii) 7f = forl = l,...,L 

iii) 7;^ = 7j^'" for I = 1, . . . ,L and z G C/R 

iv) 7f = 7f_:+ /or / = L + 1, . . . , 2L and z e C/R 

v) 7f = -7lL-i+i forl = l,...,2L 

vi) 7^^''*' = — 72L;^i for I = 1, . . . , L and z E C/M 
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Proof. Items (i) and (ii) follow immediately from Lemma 5.2 of [KSj . The other items can be 
proved as in Lemma 5.3 of |KSj if one, moreover, uses the identity T^(x, cu)"^ = J*T^{x,uj)J 
following from Wronskian identity ([7]) and invokes Corollary [1] to show that M^(ci;) is uniformly 
bounded in uo for every fixed z. □ 



5 Kotani theory 

Kotani theory links the absolutely continuous spectrum of stochastic quasi-one-dimensional op- 
erators to the set of energies with vanishing Lyapunov exponents, by using analyticity arguments 
based on a few crucial identities. In all this section it is not needed that the stochastic Dirac oper- 
ator has time- reversal symmetry or is of the particular random form given in (|30|) . Kotani theory 
for stochastic Dirac operators with bounded potentials was developed in [Sunj by providing the 
relevant identities and then following closely the arguments of [KSJ . As already mentioned, the 
paper by Sun has some obvious errors which are corrected below. Moreover, we extend the the- 
ory in order to include singular potentials and potentials which may be complex-valued matrices. 
The singular potentials model a discrete version of Dirac operators (a satisfactory discrete analog 
acting on C^'^) does not exist). 



Theorem 4 Let be given a stochastic family of Dirac operators with integrahle and singular 
potentials. Then, for k = 1, . . . ,L, the disjoint sets 

Sk = {E G M I exactly 2k Lyapunov exponents vanish at E } 

are an essential support of the absolutely continuous spectrum of multiplicity 2k. 

Just as the crucial identities are different for discrete and continuous Schrodinger operators 
(compare [KB] ), there are some variations in the formulas in [Sun] for stochastic Dirac operators 
with singular potentials as well. We need to introduce further notations in order to state them. 
Averaging over u w.r.t. P is denoted by E. Another average along the orbit of singular points is 

\ y&n[o,x] ) \ i=i / 

namely one first averages over the random sites of the singular potential. Note that Eg(l) = ps 
and that the average E can be dropped P-almost surely. Furthermore, if x§ G § is the point in 
closest to the origin, then T_xgU has a singular point at the origin and Es(/) = J P{duj) f{T_x^u). 
Hence E§ is closely linked to the Palm measure. Further the sum of the Lyapunov exponents is 
denoted by Y = ^t=i1i introduce two functions on C/M by 

= -Eg ln(det(D-M;5)) - ETt{R + M^iQ - z)) , 

and 



w 



Eg In {det{D* + B*M1)) - ETt{-R* + M'_{Q - z)) . 
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By Corollary [T] the imaginary part of is uniformly bounded away from so that w'^ are well- 
defined. The branch of the logarithm is chosen in a continuous way in z (for each uj separately) 
so that Theorem [2] then shows that is analytic. The choice of the branch is of no importance 
below. Finally for any smooth function / on we define df{uj) = dxf{T_xUj)\x=Q if ^ §. 

Theorem 5 Let '^m(z) ^ 0. 

(i) There is a constant c G M such that = wt + ic 

(ii) Y = - ^e{w'^) 

(iii) d^w' = E Tr(G^) 

(iv) 2Y = 3m(2)ETr((l + |M|n ($>m(M|))-i) 

Items (ii) and (iii) combined provide a Thouless formula for stochastic Dirac operators. The 
proof is based on a series of algebraic identities which we check first. 

Lemma 3 Let Qmlz) ^ 0. Away from singular points, the following identities hold. 

(i) aTr(ln(M^ + Mf)) = Tr(i?* + R+{Q- z){M^ - M^)) 

(ii) (9Tr((M^ + Ml)~\d^M^ - d^Ml)) = 2Tr(Gg) + d^Ti{{Q - z){Ml + Mi)) 

(iii) ±(9Tr(ln($5m(M|))) = 2 5Re (Tr(iy|)) - 53m(2)Tr((l + |M||2) (53m(M|))-i) 
where Wl = R + {Q - z)M^ and Wl = -R* + {Q - z)Ml 

(iv) d^[a^^{x,ujy'^m{Ml{T_^uj))a'^{x,uj)] = ^ $5m(z) a^(a;, cj)* (l + \Ml{T^^uj)\'^)a'^{x,uj) 

Proof. In the formulas below all functions have the argument T_^uj, and one may then set x = 0. 
Using Lemma [2](viii), a short calculation shows 

d [Ml + Mf ) = [Ml + M!) [R* - (g - z)Mi) + {R + Ml{Q- z)) [Ml + M!) . (24) 

Multiplying this by (M^ + Mi_)~^ and then using the cyclicity of the trace shows the formula of 
(i). For (ii), let us take the derivative dz of the Ricatti equation of Lemma [2]^viii) : 

d {Ml - Mf ) = - (2 + {Mlf + (Mf )2) + {Ml - M1)R* + R{MI - M't) 

+ Ml{Q - z)Ml + MKQ- z)Ml + Mi{Q - z)Mi + Mi{Q - z)Mi . 

Using this and fl24l) . some algebra directly leads to (ii) if one also uses the identity 

Tr(G^) = Tr([(M;)-i + [M^y^]-^ - {Ml + M'_)-^) , 

following from Propositions [3] and [51 

Next we turn to the proof of (iii). Let us set M|. = + lY^. with = 53m(M^). From 
M| = (M|)* follows X| = X| = (X|)* and F| = -Yl = (YD*. Straightforward calculation 
then shows 

dYl = RYl + YIR* ± Xl{Q - lke{z))Yl ± F^(Q - lke{z))Xl t '^m{z){l + {X^f - {Ylf) . 
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Thus 

dTT{\n{Yl)) = Ti{R + R* ±2X'^{Q-^e{z))) t ^m{z)TT {{Yiy\l + {X^Y - {Ylf)) 

= ±2 3fJe (Tr(iy^)) T '^rn{z)TT {{Yl)-\l + {X^Y - (Ylf) + 2Yl) 

= ±2^e {TT{Wl)) T '^m{z)Tr {{Y^)-\l + \Xl+tYl\'')) , 

where in the last step we used Tr(y^^[X, Y]) = 0. Finally let us consider (iv). When calculating 
the derivative on the l.h.s. the product rule leads to three terms. The term containing dY^ is 
given by the above formula, those involving derivatives of a^{x,uj) by Lemma[2]^vii). Hence it is 
sufficient to check 

{-R* T {Q - z)MiyYl + dYl + Yl{-R* t{Q - z)Ml) = T'^rn{z) (l + |M^p) . 
Again this follows from some algebra. □ 

Proof of Theorem (i) Set P = E Tt{R* + R + {Q - z){M^ - M^)) . By the ergodic theorem 
and Lemma [3]^i), P-almost surely 

P = lim - r dx Tt{R*{T,uj) + R{T,lo) + {Q{T,mj) - z){M^{T^uj) - Ml{T,uj))) 

1 r° 

= lim - / dx d^TT(\n{MUT_^uj) + M^{T_^u))) 
= lim - V [ln(det(M;(T_,^)+Mf(T_,^)))r^- , + 27 



y 



im rij 



where S = {xj)j(zz with xj-i < Xj and rij G Z denotes the number of branches of the logarithm 
needed in the integral from xj-i to Xj minus 1. Now by Lemma [2](iv), M^(T_(j. ._j+)Cc;) = 
M|(T_^^._^w). On the other hand, we calculate M^(T_^.^,+w) + Mf (T_^,^.+cu) by Corollary EJ^iii). 
Thus regrouping the terms shows that 

F = lim - J2 [-ln(det(Dj- -M;(T_^.^cu)Bj) -ln(det(D* + S;Mf(T_^.^xj))) + 2mnj] . 
y y 

-y<xj<o 

Hence if c is the average of 27m j over S, we have shown 

= -E§ ln(det(D - M^B)) - Eg ln(det(D* + 5*Mf )) + tc, 

and thus (i). For (ii) let us start from a formula for 7^ which follows from the identities stated 
in Proposition [8j 

7^ = lim - ln(|det(al(y,c^))|) , 

y^oo y 

where the convergence holds P-almost surely. Telescoping and regrouping gives 
Y = lim - ^ [\n{\det{a'_{xj+,u))\) -\n{\det{a'_{xj_i+,uj))\)] 



lim — 



ln(| det(a!l(xj+, cij)a^(a;j — , cu) "'^)|) + / fix ln(| det(Q;l(x, u; 



18 



The first contribution can be evaluated with Lemma [2](vi) and the definition of Eg, the second 
contribution be summed up and the integrand evaluated: 

Y = -E§ ln(|det(D* + B*M!)|) + JJe lim - F dxi:i{a'_{x,oj)-^d^a'_{x,u)) . 

y^oo y Jo ^ 

Finally the last expression can be calculated using Lemma Et^vii) and then the ergodic theorem 
completes the proof of (ii). Let us point out that one could have started from 

Y = - lim -\n{\det{aliy,uj))\) . 

y^oo y 

Then a similar calculation leads to 7^ = —^e{w^). Because are analytic, this also provides 
an alternative proof of (i). 

(iii) Let us set J' = 2ETT{G^)+d,ETT{R-R* + {Q-z){M^+M^)). Because the probability 

of having a singular potential at vanishes, E Tr(G^) can be replaced by E Tr (Gq) . Furthermore 
the term R — R* drops out due to the derivative dz- Hence Lemma [3](ii), the ergodic theorem 
and reordering of the terms imply as above that P-almost surely 

r= hm - V Tr((M;(T_,a;) + Mf(T„,.cu))~^(M^(T_,,cu)-M!(r„,.cu)))|^^- , 

where we also used the left-continuity of a; G M h-* M^(Tj.cc;). The terms with Xj— now have 
to be evaluated using Lemma [2](iv) or its equivalent formulations. The factor (M\(T_x.j^iS) + 
M^_{T_x^j^Lo))^^ is given by the inverse of Corollary [2]^iii). Corollary [2]^ viii) moreover allows to 
calculate M^(T_a;^.+ti;) — Mf (T_i.^.+c<j) . Replacing both identities then shows 

.r = EsTt (^{Ml + Ml)-^ {D + M'_B){D ~ M^B)-^ Ml 

+ Ml {D* + B* M1)-\B* Ml - D*) - + ) 
= dzEs [Tr(ln(D* + 5*Mf)) - Ti{\n{D ~ M^B))] . 

Due to the definitions of and this concludes the proof of (iii). 

(iv) We set = E(2 3fJe(Tr(iy|)) - ^m{z)Tr{{l + |M||2)S5m(M|)-i)) . By Lemma EKiii) 
and the ergodic theorem one has P-almost surely 

±Kl = lim - y ln(det{Qm{Ml{T_^Lj))))p~ . 

Now evaluate 3m(M^(T_^^.+u;)) by Corollary [2](vi). This implies 

Kl = ^E§ ln(det(|DTM^5n) = ^ 2 3fJe E§ ln(det(D ^ M|5)) . 

Similarly, using Corollary El^vii), = ^2 3fJeE§ ln(det(D* T B*Ml)). From these identities 
one readily completes the proof. □ 



The second part of the following theorem establishes Theorem 6.6 of |KSj also for complex 
valued potentials. 
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Theorem 6 Consider the positive operator f/|. = (5>m(M|))2(l + |M|n-i(5>m(M|))2 and 
denote its eigenvalues by > ... > u\^± > 0. Further let E E M., e > and k = 1, . . . ,L. 
Then 



1 9 



(25) 



If furthermore E is such that 7/^ = liiiie^o 7;^'*'*'^ exists for I = 1, . . . , L, then 



k 



1=1 



< 



1=1 



=k+l 



(26) 



Proof. This is an adaption and shght generahzation of the proof of Theorems 6.5 and 6.6 of 
[KSj (the reasoning in |Sunj is erroneous at several points). For any L x L matrix F let A^F and 



dA F the second quantizations on the fermionic tensor product AC, such that e 
Let z = E + ie and F| = 3m(M|). Define Fl{x,uj) = F|(T_^w)^a^(x, w). Then 



A^e 



d^A^\Fl{x)\^ = A'^Flix)* {dA''{Fl 



d,\Flix)\' Flix)-')A'F.^' 



X) 



Thus by Lemma [3](iv) 

d^A''\Fl{x,uj)\^ = T'^rn{z) A^Fl{x,ujy {dA^Ul{T_^u})-^) A^F__ 
so that for ^miz) > 



d^A''\F^{x,uj)\'' > - Qm{z) WdA^U^iT^^uy^W A''\F" 



d^A^lFlix.uj)]' < "^miz) WdA^UliT^^iuyW A''\Fl{x,uj)\' . 
Integrating hence gives 

A^|F^(x,cj)|2 > exp(^-Qm{z) dy\\dA''Ul{T_yUj)-^\\^ A''\Fl{0,uj)W 
A^|Ff(x,cj)|2 < exp(^Qm{z) dy\\dA''Ul{T_yUj)-^\\^ A''\Fl{0,uj)\^ . 

Note that by Lemma[2](vi) and Corollary [2]^ vii) the functions \F^{x,uj)\ are actually smooth also 
for X G S. We combine this with the inequalities 



\\A''aUx,u)f\\A'Yl{T^,LU 
Taking logarithms thus shows 
||AV(x,cu)f 



ri||-i < ||A'=|F^(x,cu)n 



< ||A'=a^(x,o;)||^ 



lA'^YliT.^u) 



In 



\A>'Yl{T_^uj)-^ 



< ismiz 



r dy\\dA''U^{T.yu;)-'\\ + In (||A'=|F^(0, a;)p 
Jo 



Now by Corollary [H Y^{uj) ^ is uniformly bounded in u. Thus dividing by x and then taking 
the limit x 00 shows by Proposition [Ht^iii) and the ergodic theorem 



2 ^7f < ^m{z) E \\dA''{U_ 



ktTTZ\-l\ 



'^m(z' 



^ 1 

E y - 

l=L-k+l - 



(27) 
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Combining this with k replaced hj L — k together with Theorem [5](iv) stating 

L L ^ 

1=1 1=1 
proves inequahty fl25|) for the sign — . Similarly one has 



ln(||A'=a;(x,cj)f ) > - !3m(;2) [ dy\\dk^Ul{T_yUJ 

Jo 



II ^ \\\A^Y^{T^^uj)\\ 



As the last term is bounded along the orbit, Proposition [8](iii) now implies 



1=1 l=L-k+l 

which again combined with Theorem [5]^ iv) proves fl25|) for the sign +. 

For the proof of fl26|) we need the following general fact. If T, 5* > are two positive ma- 
trices, then the positive operators T2ST2 and S^TS^ have the same spectrum (this follows 
from Tr{{T2ST2)^) = Tr((^2r^2 )'^) for all n G N). Hence are also the eigenvalues of the 

imaginary part of the Herglotz function (1 + |M|p)~^M^(l + |M|.p)^^ and by the Herglotz 
representation theorem it follows as in [KSj that 

^ > -^-^ for e > 5 > . 



Combining this fact with Theorem [5](iv) and the bounds ( 1271) and (l28ll gives 

1=1 "±,/ 1=1 "±,Z i=k+l ^±,1 1=1 1=1 

Now taking the limit 5 — leads to (126|) . □ 
From this point on the proof of Theorem H] is line by line the same as in (KSj . 



6 Time reversal symmetry and Coupling Hypothesis 

None of the results of Sections [2] to E] used the time reversal invariance ([2]). In this section, we 
first implement this symmetry and then describe the model of Theorem [1] in more detail and 
state the Coupling Hypothesis. The proof of the following result is immediate. 

Proposition 9 Suppose that H is time-reversal invariant, namely satisfies ([2]). Then 

TY^J = T^{x) , Wl = -(Ml)-' . 
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Proposition 10 If has time-reversal symmetry, the averaged Green matrix satisfies 



TG\,J = (G^)* , J*^m{G'^)J = 53m(G^) . (29) 

If furthermore (j) = {v, Jv) for some v G C^^ satisfying v*Jv = 0, then the 2x2 matrix (f)*G\;(j) 
is a constant multiple of the identity. 

Proof. The Hamiltonian satisfies J*HvJ = ^ so that J*{Hv - z)-^J = {Hv-z)-^. This 
implies that for any vectors v,w E C^-^, v* J*G^Jw = W*G^v = v*{G^yw which implies the 
first identity in fl29l) . from which the second one can be directly deduced. As to the last point, 
for any vector w one has w*G^w = w*{GyYw = w^J*G^Jw = {Jw)*G^Jw. Moreover, for 
any w = \v + X'Jv G Ran(0), one checks the orthogonality w*Jw = 0. These facts imply 
w*G^w = i Tr(0*G^0) llwf . □ 

The last statement of Proposition [10] reflects Kramers' degeneracy stating that the spectrum 
of a time reversal invariant Hamiltonian with odd spin has even multiplicity. In particular, 
for eigenstates Hip = Eip gives HJ'ip = EJ'ip. For the same reason, the singular values of the 
transfer matrices are degenerate (see Lemma ID^ii)) which implies the degeneracy of the Lyapunov 
spectrum. 

Next let us come to the construction of the stochastic Dirac operators of Theorem [T] and of 
the associated dynamical system. Let s G [0, 1) = M/Z. Each operator II{uj) is of the form ([T]) 
with singular potentials at S = Z + s, hence xj = j + s. The Vj are drawn independently and 
identically out of J'so*{2L) with some probability law pv with compact support. Furthermore 
the potential W G Jso*{2L)) is of the form 

K 

^i^) = "^'^(^ + S-J + 1), (30) 

jez k=i 

where K eN, each Wk G L^S^, Jso*{2L)) has support [0, 1] and the vectors {Xj,k)k=i,...,K e 
are also drawn independently and identically according to a probability distribution pw with 
compact support. Then is a compact subset of {J'so*{2L) x R^)^^ x M/Z and P = (pw x 
Pv)^^ X ds. The R-action T is the natural right shift on Q and P is indeed ergodic and even 
mixing w.r.t. T. In order to state the main hypothesis on the randomness, it is convenient to 
introduce the transfer matrix T^(W, V) as the solution T^(l, 0) of (jS]) with potential W and jump 
e-^^ at 1. Setting Xj = {Xj^k)k=i,...,K (which determines the potential bump Wj = ^A^i-^i-fc^fc 
between j — 1 and j), this notation implies T^{Xj, Vj) = T^{j + s, j ' + s — 1, tu) where the transfer 
matrix on the r.h.s. is defined by ([5]) with the Hamiltonian II{uj). 

Coupling Hypothesis: The semi-group generated by {T^{X,V) \ {X,J'V) G supp(p>v x Pv) } 
is Zariski dense in S0*(2L) for all G R. 

Let us stress that this hypothesis can be verified if pw x Pv is supported on a finite set of 
points, and also if either pw or pv is concentrated on a single point, notably the disorder is given 
only by a random potential W or the random Dirac peaks Vj6j. Furthermore this hypothesis is 
satisfied whenever the set of T^(A,V) contains an open set (this property does not depend on 
E). This is e.g. the case if pv contains an absolutely continuous part w.r.t. to the Haar measure. 
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7 The Lyapunov spectrum 



This section proves a criterion for the distinctness (apart from Kramers' degeneracy) of the 
Lyapunov exponents for random products of matrices in S0*(2L). It can be immediately apphed 
to the transfer matrices if the Couphng Hypothesis holds. On the other hand, we believe it to 
be of somewhat independent interest and thus took care to make it readable without reference 
to the rest of the paper. Instead of the group S0*(2L) as defined in the introduction it will be 
more convenient to work with an isomorphic group G for which the polar decomposition takes a 
more simple form. Thus we define in case of even L = 2d and odd L = 2d + 1 respectively 



where d x d square matrices carry the index d. Then introduce A = dmg{A, A) which satisfies 
A* = A~^ and set G = ^*S0*(2L)^. This group consists of all 2L x 2L matrices M satisfying 



where S = dia:g{A^A,A^A). Note that the matrices and S commute, J^* = — J = J ^ and 

5* = 5 = s-^. 

Lemma 4 Let M ^G and v G C^^. 

(i) M* eG 

(ii) IfMv = \v, then M*Jv = \-^Jv, MJSv = XJSv and M*Sv = T^Sv. 

(iii) The vectors v and JSv are linearly independent for f 7^ 0. 

(iv) For M> e G, there exists U eGf] SU(2L) such that UMW = V, where 
V = diag(ai, . . . , a^, 1, a^^, . . . , a]^^, . . . , a^-^, 1, oi, . . . , a^) if L = 2d + 1 and 

T> = diag(ai, . . . , a^, a^^, . . . , a^^, a^^, . . . , a^^, ai, . . . , a^) in case L = 2d, with real constants 
0-1 > > • • • c^d > 1- Note that V & G. 

(v) There are unitary matrices /C, W G G fl SU(2L) and a diagonal matrix T> as in (iv) such that 
M = fCVU. 

(vi) One has det(A^) = 1 and the group G is connected. 

Proof, (i) follows by inverting the relations in (13T!) . For (ii) note that M*JM. = J implies 
J*M*J = M~^. Hence J*M*Jv = X'^v implies M*Jv = X'^v. From M^SM= S it follows 
that SM^S = M'^ = J*M*J. Taking the transpose one obtains SMS = J*MJ and hence 
SM.SJv = —XJ*v and therefore M.SJv = XSJv. Now using the same calculation as above 
yields the last equation. 

(iii) Writing v = (^) and J'Sv = Xv gives Xa = A^Ab and Xb = —A^Ad. As A* A is real and 
[A^AY = 1, this implies |Apa = A^AXb = —a and therefore (1 + |Ap)a = implying a = and 
6 = and hence v = 0. Therefore these vectors are linearly dependent if and only if f = 0. 

(iv) First we need some basic facts. We say that a subspace V of C^^ is G-like if for any vector 
f G V one has J'v,Sv, J'Sv G V. The space spanned by v,J'v,Sv and J'Sv is G-like. The 




M*JM 



J 



M'SM 



(31) 
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intersection of two G-like subspaces is G-like. Furthermore, if V is G-likc, then also the orthogonal 
complement V""- is G-like. To see this, take v eY,w E V"*- then {Sw, v) = (W, Sv) = {w, Sv) = 0, 
and {Jw, v) — —{w, Jv) — 0. Therefore Sw^ Jw e V"*- and hence also JSw e V"*-. 
For > the eigenspaces are orthogonal. Let Vi be the cigcnspace for the value 1 (possibly 
only the zero vector) and Vq be the orthogonal complement. By (ii) and the consideration above, 
these spaces arc G-like and they are invariant under }A*}A. By (ii) and (iii) the dimension of 
Vo is divisible by 4, say dimVo = 4r. 

First claim: Vq has an orthonormal basis of eigenvectors of M*M. of the form 

Vi, V2,..., Vr, Slh, Sv;, JVi, JVr, JSvi, JSv^. 

Indeed, if dim(Vo) = 0, there is nothing to prove. Otherwise let > 1 be the biggest 
eigenvalue of Ai*Ai which is also the biggest eigenvalue of Ai*Ai restricted to Vq and let v e Vq 
be some corresponding eigenvector. Then J'Sv is another eigenvector for the same eigenvalue. 
Take w = v + fiJSv, where G C can be chosen in such a way that w and JSw are orthogonal. 
Then also Jw and Sw which are eigenvectors to the eigenvalue a^^ are orthogonal. As Oi > a]~^, 
the space spanned by w and JSw is orthogonal to the space spanned by JW and Sw. Therefore 
normalizing w to Vi = w l\\w\\ the vectors Vi, Svl, Jv\, JSvl are orthonormal. Denote the space 
spanned by these vectors by Vo,i C Vq and its orthogonal complement in Vo by Vo,2 which is 
again a G-like, Al-invariant subspace. One proceeds by induction to complete the proof of the 
claim. 

Second claim: If L = 2d, then dim(Vi) is divisible by 4 and there is an orthonormal basis 
of the form Vr+i, . . . ,Vd, Su^^, . . . , Sv^, Jvr+i, ■ ■ ■ , Jvd, JSv^^, . . . , JSvd- If L = 2d + 1, then 
dim(Vi) is congruent to 2 mod 4 and one has an orthonormal basis which is of the form 

Vr+l, ■ ■■,Vd, Vd+l, 5tV+r, • • • , Sv2, JVr+l, ■ ■ ■ , Jv,i, JVd+l, JSu;:^, . . . , JSv2 with SVd+l = Vd+l- 

Indeed, as J is unitary and operates on Vi, there is an orthonormal basis of Vi of eigenvectors 
of J . The eigenvalues of J are If Jv = ±zv, then JSv = SJv = ^zSv. Hence the 
dimensions of the eigenspaces of J in Vi are equal. If dim(Vi) > 4, there are two orthonormal 
vectors wi,W2 satisfying Jwj = iWj. As JSw] — —iSw] the vectors Wi,W2,Swi,Sw2 are 
orthonormal. Set v^+i = "^(''^i + SW2). Then the vectors Vr+i, Jvr+i — -^{wi — SW2), Sv^^ = 
■^{^2 + Swi) and JSvr+i = :^('^2 — Swi) are orthonormal. They span a 4-dimensional G-like 
subspace of Vi. Denote its orthonormal complement in Vi by V2 and proceed by induction to 
obtain the vectors Vr+2, ■ ■ ■ ,Vd. In case L — 2d this shows the above claim; if L — 2d+ 1, one is 
left with some 2-dimensional, G-like subspace Yd-r+i- This space is spanned by the orthonormal 
vectors w and Sw where Jw — iw. Set Vd+i — ;^('"^ + then Vd-^i and Jvd-^i form an 
orthonormal basis of Nd-r+i and Svd+i = Wd+i. 

Construction of U: From the first two steps we obtain an orthonormal basis of eigenvectors of 
M*M of the form {vi, . . . ,Vd, Vd+i, SW, ■ ■ ■ , SvJ, Jvi, . . . , Jvd+i, JSvl, • • • , JSvt) in case L = 
2d-\-l and the same without the entries containing Vd+i ii L — 2d. The corresponding eigenvalues 
ol vi, . . . ,Vd shall be denoted by of > a| > . . . > 1. The eigenvalue corresponding to Vd+i if 
L = 2d -I- 1 is 1. Denote the canonical basis of C^^ by e^, i — 1,. . ., 2L. Let us define the unitary 
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matrix lA by 



L = 2d L = 2d + 1 

Uvi = Cj i = 1, . . . ,d Uvi = Cj i = 1, . . . ,d + 1 

USWi = ei+d i = l,...,d USWi = Ci+d+i i = l,...,d 

Ujvi = -ei+2d i = l,...,d Ujvi = -ej+2d+i i = l,...,d+l 

lAJSWi = -Cj+sd i = l,...,d UJSvl = -ei+3d+2 i = l,...,d . 



Then defining the diagonal matrix V as in the statement of the proposition, one has lAMlA* = V. 
For i = 1, . . . ,d, one has 

{WJU)vi = WJe, = -Wci+L = Jvi , 

{WJU)Jv, = -WJe,+L = -U*e, = -v, = J{Jvi) , 

similar calculations hold for Sui, JSvl and also Vd+i, Jvd+i in the case L = 2d + 1. Thus 
one obtains U*JU = J . It is a matter of calculation to verify that U^SU = S and hence 
e G n U(2L). Finally, as W G G we have AUA* G S0*(2L) n U(2L) = SP(2L, M) n 0(2L) and 
hence det(W) = det(^^*) = 1 and therefore U G SU(2L). 

(v) As M*M G G and M*M > 0, by (iv) we find W G G n SU(2L) and a diagonal matrix 
V as above, such that UM*MU* = V^. Set /C = MWV'^ G G, then M = KVU and 
/C*/C = V-^UM*MU*V-^ = 1. Hence JCeGH U(2L) = G n SU(2L). 

(vi) By (v), det(A^) = det(/C) det(P) det(W) = 1. Furthermore as the group SP(2L,M) n 0(2L) 
is connected, also S0*(2L) fl U(2L) is. Using the decomposition in (iv) one easily obtains that 
G is connected. □ 

Now let (3^n)n>i be an i.i.d. sequence in G. Then by Lemma H] the whole associated Lya- 
punov spectrum has at least multiplicity two. So let 7i, 7i, 72, 72, • • • , 7l, 7l be the 2L Lya- 
punov exponents with 7i > 72 > ■ ■ • 7l- Lemma H] also shows 7p = —•jl+i-p and in the case 
L = 2d + 1, one has 7^+1 = 0. Therefore it is always enough to consider 71,..., 7^. Set 
v^^ = ei A . . . A Cp A e2L^d+i A . . . e2L-d+p and define Lp = spanuHA"^^ Mv^-^^ \ M G G}) which 
is a real linear subspace of A^^C^^. Note that hp does not have to be a complex vector space. 
Taking the real part of the scalar product on A^^C^'^ induces a scalar product on Lp but actually 
one does not need to take the real part as the following lemma shows. 

Lemma 5 The scalar product in A^^C^^ of two vectors in hp is real. Let fi, /2, /s, /4 G hp and 
consider /i A/2, fz^fi on one hand as elements in A^(A^^C^'^) and on the other hand as elements 
in A|Lp considered as tensor product over the field M. Then the scalar products coincide, i.e. 

(/l A /2, /s A /4)A2(A2pC2i) = (/l A /2, /s A /4)a2Lp- 

Proof. One finds JSci = —e2L-d+i and JSe2L-d+i = e, for i = 1, . . . , c? which implies 
K^P{JS)v^P^ = = v^Pl For M eG one has SMS = J*MJ and hence 

{v^p\A^PMv^p'^) = {k^'PSv^P\k^-p{SMS'^)v^P^) 

= {A'^f{JS)v^P\A'^P(MJS)v^P^) = {v(p\A^pMv(p')) . 
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Therefore {K^pMv^p\ A^^Mv^^^) = {v'-p\ A^^{M*M)v^p^) is real for all A^,7V G G and by linearity 
the h?P£?^ scalar product for two vectors in Lp is real. The second statement follows from the 
first one using {f\ A /a, h A U) = {fi, h){f2. U) - (/i, A) (A, h). □ 

Considering fi A /2 as element in A^Lp on one hand and as an element of A^(A^^C^^) on 
the other hand induces an M-linear map A^Lp — > A^(A^^C^'^). By Lemma [5] this map preserves 
the inner product and is hence injective. Therefore A|Lp can be viewed as real subspace of 
A^(A^^C^^). The following criterion for distinctness of the Lyapunov exponents is adapted from 
[GHllBL] . 

Definition 1 A subset T of G is hp-strongly irreducible if there does not exist a finite union W 
of proper linear subspaces of hp such that (A^^A^)(W) = W for any M. in T. 

Proposition 11 Let {yn)n>i be a sequence of i.i.d. random matrices in G for L = 2d or 
L = 2d + 1 and let p be an integer I < p < d. Let T be the semi-group generated by the support 
ofyn- Suppose that T is 2p- contracting and hp-strongly irreducible and that E(log_,_ ||3^i||) < oo. 
Then jp > 7p+i. 



Proof. Let k be the dimension of L, 
on. For any A4 G 



^ ,p and 

let A4 denote the matrix in G1(A;, 



an orthonormal basis to be chosen later 
R) with the entries 



M 



{f,,A'^Mf,), l<t,j<k 



If W G G n U(2L), then A^^W G A^^G n ViA^PC"^^) and hence the restriction of A^pU to Lp is 
orthogonal, i.e. U E 0(Lp). Let us use the notation A'^^M = ^{M) . One has \\M\\ < \\^{M)\\ 
as Lp is a subspace of A^^C^-^ and by Lemma [5] one also obtains ||A^A<|| < ||A^$(A^)||. 

> flrf > 1 be the singular values of Ai as occurring in the de- 



and ||A2$(;W)|| > \\AIM\\ > 



Claim: Let Oi > 02 > 

composition in Lemma Hl^v), then = a1---ap 

■ a\ - ■ ■ ctp_iCtp+i. In the case p = d, we define a^+i = a^^. 

Indeed, set /i = v^p'> = ei A . . . A Cp A e2L-d+i A ... A e2L-d+p and if p < ci set /2 = ei A . . . A 
ep_i A Cp+i A e2L-d+i A ... A e2L~d+p~i A e2L-d+p+i- In the case p = rf, set /2 = ei A . . . A e^-i A 
CL+d A e2L~d+i A ... A e2L-i A e^. Further, for any d x d invertible matrix B and any matrix C 
with B*C = C*B, one can construct the following element of G: 



AT 






c 







cos(v9) 



- sin((y9) 



















{B* 



,-1 










sm{(p) 



cos{ip) 




\ 



-c 




B / 



if L = 2d pencil out 
the rows and columns 
containing if. 



(32) 



Thus for p < d, one readily finds A/" G G with /2 = A^pM fi G Lp. In the case p = d define A^i by 
setting B = 1 and Cij = except C^^ = 1 and define A/2 by setting 5 = 2-1, C = 0. Then one 
obtains (22('^-i) - 22(^-2))/2 = (22(^^1) - A'^PM2){A^PMifi - /i) G L^. In conclusion, /i,/2 G Lp 
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can be completed to an orthonormal basis of hp. Now let us write A4 = ICDU as in Lemma Hl^v) , 
then 



||$(7W)|| = al---al = ||A'Pp/i|| = < \\V\\ < \\A^PV\\ = \\<^{M)\\ 

where the last inequality holds as Lp is a subspace of A^^C^^. Hence = but 

\\V\\ = \\ICVh{\\ = \\M\\. As mentioned above, ||A^$(A<)|| > ||A^A^||. Furthermore one has 

\\A'M\\ = \\A'V\\ > ||A2%iA/2)|| = \\M\\-al---al_,al^, 

Hence the claim is proved. 

Let T be the semi group induced by the distribution of As T is Lp-strongly irreducible, 
clearly T is a strongly irreducible subset of G1(A;,M). As T is also 2p contracting, there exists a 
sequence (A^„)„>i in T such that lim„^oo ||'^'(A1n)in|A^$(A^„)||~^ = oo. As \\Mn\\ = ||$(A^„)|| 
and ||A^$(A^„,)|| > ||A^A^„|| by the above claim, one obtains 

lim \\M4^\\A^MnV > lim \mMnW\\A^HMn)V = oo. 

n— >oo 71— +00 

Hence T is contracting. The two biggest Lyapunov exponents associated to the sequence (3^ri)„>i 
shall be denoted by 71 and 72. Then by the claim, the definition of Lyapunov exponents and 
(BLj A.HI.6.1] one has 

p p-i 
2 J]^^ = 71 > 72 > 2^7i + 27p+i , 

1=1 2=1 

implying 7p > 7p+i. By definition of Op+i one actually would have to replace 7p+i by 7p+2 = 7d+2 
in the case L = 2d + l,p = d. Then one gets ■jd > ld+2 = —Id and therefore 7^ > = 7d+i- n 



Theorem 7 Let (3^n)n>i &e a sequence of i.i.d. random matrices in G for L = 2d or L = 2d + 1 
Let T be the semi-group induced by the support o/3^i and let E(log_|_ ||3^i||) < 00. Suppose that T 
is Zariski dense in G, then all Lyapunov exponents are distinct. 

Proof. According to the proof of Proposition [TT] the inequality 7p > 7p+i follows from the fact 
that the semi-group T = {Ai | 7V1 G T} is strongly irreducible and contracting in GL(/c,M) as 
defined above. Now T is Zariski dense in G = {Ai \Ai G G}. Otherwise there would be a 
polynomial P on GL(fc, M) such that P(T) = and P{Ai) 7^ for some G G. As the entries 
in A4 are polynomials of the entries in A4, this leads to a polynomial P on GL(2P, C) such that 
P(T) = and P{A4) 7^ for some Ai E G, contradicting the fact that T is Zariski dense in G. 

Now suppose T is not strongly irreducible. Then there would be a finite union of proper 
subspaces W = Vi U . . . U V„ such that A^(W) C W for all M eT. The property M{Yi) C 
can be written as {w,Aiv) = for all w G V^,f G Hence the set of all such matrices At 
is Zariski closed. The property At{W) C W is therefore a finite intersection of finite unions 
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of Zariski closed sets and hence Zariski closed. As T is Zariski dense in G, this then implies 
G(W) C W. Therefore, if G is strongly irreducible, then also T is. 

To show that T is contracting we want to use Theorem 6.3 of |GMj which states that if the 
algebraic closure of T is strongly irreducible and contracting, then also T is contracting. Hence 
it is only left to show that G is strongly irreducible and contracting. 

The property of G to be strongly irreducible is equivalent to G being Lp-strongly irre- 
ducible. As G is connected we have to show that there is no proper subspace V C Lp such 
that {A'^PM){Y) C V for all M e G. Suppose such a V exists. For ai > 02 > . . . > > 1 take 
V = diag(ai, . . . , a^, 1, . . . , \ a]^^, . . . , a^^, 1, ai, . . . , a^). The relation (A^^P")(V) C V im- 
plies that either v^^'^ G V, but then Lp = V or that v^^ is in the orthogonal complement V"*". But 
then by Lemma H^i) one has, for v e Y and any M e G, {A^pMv^p\v) = {v^p\A^pM*v) = 0. 
Hence Lp = V"*". Therefore V is not proper. 

Now it is only left to show that G is contracting. By the proof of Proposition [11] this 
follows if G is 2p-contracting. Therefore take a matrix Ai of the form fl32l) with C = and 
B = diag(Ai, . . . , A^). such that all moduli of the eigenvalues are distinct except for the fact 
that always two eigenvalues have the same modulus. The sequence AI" then shows that G is 
2p-contracting. □. 

Proof of Theorem [T](i) and first claim of (ii). The Coupling Hypothesis implies by The- 
orem [7] that the Lyapunov exponents as defined in Section [6] are distinct apart from Kramers' 
degeneracy. The symplectic symmetry of the Lyapunov spectrum implies that no Lyapunov ex- 
ponent vanishes for even L, while for odd L there are exactly two vanishing Lyapunov exponents. 
By Theorem [5] the absolutely continuous spectrum is absent for even L and has multiplicity 2 
for odd L. □. 



8 Absence of singular spectrum 

In this section we only consider the random model described at the end of Section HI For any 
configuration u = {{\j,k)k=i,...,K;j&, iyj)jez, s) E Q let Co denote u excluded the singular potential 
V = Vo at s, i.e. o) = {{\j±)k=i,...,K,j£Z , (VjOjgZjyo, s). The distribution of cD shall be denoted by 
P and that of V by py. With these notations P = P x pv. We only consider the case where L is 
odd and pv is absolutely continuous w.r.t. to the Lebesgue measure. Next recall the definition 
f|T8l) of V G J'so*{2L). Note that V is only defined for almost every V and for almost every V 
there is a pre-image V, which is not necessarily unique. Furthermore the pre-images of zero sets 
are zero sets and hence the distribution p^ of V, i.e. the image measure of pv, is absolutely 
continuous w.r.t. the Lebesgue measure on the vector space j7'so*(2L). 

As V denotes the singular potential at xq = s, let Gy denote the averaged Green matrix 
at the point xq = s, that is, Gy = Gy{s) with the notations of Proposition [51 Note that this 
matrix actually depends on u = (u;,V), but in most of the arguments below oj will be fixed. 
Furthermore, Proposition [7] shows that Gy actually only depends on V (which is a real statement 
statement since the map V 1— >■ V is not injective). Hence it is sufficient to prove almost sure 
statements w.r.t. the distribution py of V instead of w.r.t. the distribution pv of V. 
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Let = /ij),v denote the associated positive matrix valued measure. The function E ^-^ 
is in L^{fii^) for all u. On the set of such measures one may introduce the weak-* topology induced 
by the functions E i— > Qm{{E ~ z)^^) for z in the upper half plane. As the pairing of this function 
with the measure fi^^ is just '^m{G^), it follows that the map u ^ fi^ is Borelian. Finally let 
fJ'oj,k = yUw,v,fc denote the measure corresponding to elG^Ck where is the fc-th canonical basis 
vector of C^^. 

The aim of this section is to prove that almost surely in uj the measure fi^^ is absolutely 
continuous or equivalently, that its singular part vanishes, i.e. fiuj,sing(M) = 0. Therefore we will 
first show that almost surely one only needs to consider /x^^^i and then we show that /i(^,i,sing(IR) = 
almost surely. To obtain the first part we compare the measures /ii>,v,i and ficj,v,k for fixed uj and 
show that they are almost surely equivalent. Once cyclicity issues are settled (Proposition [T2l) 
and matrix analogues of rank one perturbation results are proved (Proposition [T3l) . the proofs are 
basically modifications of the arguments of [JLj . Our starting point are the following observations 
linked to Kramers' degeneracy. 

Lemma 6 For 1 < k,l < L let us introduce the 2L x 2 matrix = (e^, 6^+^). 

(i) Let j denote the 2x2 symplectic form, then J^k = ^fcj - 
Furthermore one has ^^^fc e Jso*{2L) and ^fej^^ + ^;j*^fc G Jso*{2L). 

(ii) For yi, 3^2 e Jso*{2L) one has 3^i3^23^i e J'so*(2L). 

(iii) '^IG^'^k is a multiple of the unity matrix, which means '^IG^'^k = ^XG^eu 1. 

Proof. The identity J'^k = ^kj is readily verified. Furthermore (\E'fc\E'^)* = '^k'^l and one 
has J*-^k'^lJ = = ^fc^fc = (^fc^fc)* showing ^^^1. e Jso*{2L). Similar calculations 

show 'i'kj'^^ + ^d**fc e Jso*{2L) e Jso*{2L) and (i) is proved. To obtain (ii), first note that 
3^1,3^2 are self-adjoint and hence 3^i3^23^i is self-adjoint. Furthermore one has JT"*}^! 3^23^1 17 = 
J*yiJJ*y2JJ*yiJ = yiyiyi = (3^13^23^1)* and also (U) is proved, (iii) is just a special case 
of Proposition [101 □ 

The measure class of /i^ is given by the trace, i.e. by the sum Yl'k=i t^^,k = 2 J2k=i t^i^Mi where 
the last identity follows from Lemma EJ^iii). 

Proposition 12 For fixed Cj, one has that for Lebesgue almost all V G j7'so*(2L) the set of 
energies {E G M | Gy"*"*" exists and ^JG^'^k is invertible } has full Lebesgue measure. 

Proof. We first claim that for fixed z in the upper half plane Ui, there is a V G J'so*{2L) such 
that is invertible. Recall that G^ = ((Gg)"^ + V)-\ Set (Gg)-^ = X - ty-^ with 

37-1 = -^m{{G'o)-^) > 0. As J*G'oJ = (Gg)*, one has X,y~\y G Jso*{2L). Then consider 
V = -3fJe((Gg)-i) + XV with a perturbation V G Jso*{2L). Then 

G^ = (-^37-1 + AP)-i = ly + xyvy-tx^yvyvy + oix"^). 

Note that V now depends on A and V, furthermore 3^'P3^ G j7'so*(2L) as well as 3^'P3^'P3^ G 
Jso*{2L) by Lemma El For any 2 x 2 matrices A, B, G one has det(A + XB + X^G) = det{A) + 
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XTT{A{j*Bjy) + (det(fi) + Tr(A(j*Cj)*)) + 0{X^). Furthermore for W G Jso*{2L), one has 
{j*^*W^kjy = j**fcW**,j = '^lJ*W^J'^i = ^^W^/. Thus from the above 

det(^;G^^,) = t det(^;3^^fc) + tXTT{^*iy^,^iyvy^i) + 

+ A' (det(*;3^P3^^fc) - z Ti{^*iy^k^iyvyvy'^i)) + ^(a^) . (33) 

If det(^;y^fc) ^ 0, then the claim is true (just take A = 0). If det(^;3^^fc) = 0, but "^jy^k ^ 0, 
then set V = y-^ E Jso*{2L) and ([33]) reduces to 

det(^;G^^fc) = «ATr((^;3^*,)*(^,m)) + ^(A^) . 

Since the coefficient before A only vanishes if \E'^3^\I'fc = 0, this is not equal to zero for small A and 
the claim holds again. Finally, if ^^D^^fc = 0, then set V = 3^-^(^ij^^ + ^^j*^; + "ifi^Dy-^ 
which lies in J'so*{2L) by Lemma El part (i) and (ii). Then (1331) reduces to 

det(^;3^^fc) = AMet(^;^j^^^fc + ^K^kf + ^iWi^k) + O(A') = A^ + ^(A^) , 

where we used \E'^\E'a; = 5/,^- Hence this determinant is again not zero for small A. Thus for all 
cases we find some V such that ^^G^^k is invertible and the claim is proved. 

Now by definition of the determinant and Cramer's rule the function V ^— det(\E'^Gy\E'fc) = 
det(\&*(((jQ)~^ + V)~^\E'fc) is a rational function on the vector space j7'so*(2L) which does not 
vanish completely by the claim above, therefore it does not vanish for Lebesgue almost every 

V G J'so*{2L) w.r.t. the Lebesgue measure on J'so*{2L). 

Next recall that the boundary values Gy"*"*" exist almost surely in E by analyticity. For 

V as described above, the map z i— > det(\E'^G'y\I'jt) is analytic in the upper half plane and 
does not vanish identically. Therefore for Lebesgue almost every E, exists and one has 
det(^;G:^+*°^fc) ^0. □ 

Proposition 13 Let uj and V G Jso*{2L) be fixed and define Va = V + A\I'fc\E'^. 

(i) The set Av^,k = {E eR\ ^^G^+*°^fe exists and 55m(^*G^+'°^fe) > 0} is independent of X 
and it is an essential support of the absolutely continuous part of ^cb,Vx,k- 

(n) The singular part of iJcj,Vx,k is supported on the set {E G M | '^IG^^^^'^k = — A^^ 1}. 

(iii) For any i? C R of zero Lebesgue measure, we have ^JiciiVx,k{B) = for Lebesgue a.e. AG M. 

Proof, (i) We prove that A^ ^ = ^Vo,k C ^Va.A; ^ill A; the other inclusion can be obtained 
analogously. Hence let E G A^^k- We first claim that 1 + A\E'fc\I'^Gy^*° is invertible. Suppose 
(1 + X'^ k'^lG^~^^'^)v = 0. Then v is in the range of and there are a,?7 G C such that v = 
aek+(3ek+L. We use elG^^'^CL+k = = el^^G^+'^'ek following from J*G^yJ = (G^)*. Thus a = 
-XaelG^+'^Ck and (3 = -X/Sel^^G^+'^ek+L- But as 53m(e^G:^+*°efc) = ^m{el^^G^^'^ek+L) > 
for E G A^ k this implies a = = P and hence v = 0. Therefore the kernel of 1 + X^ k^l^y^^^ is 
indeed trivial. Hence by Proposition [H G^+*° = G^+*°(l + A^fc^^G^+*°)-i exists. Furthermore, 
also by Proposition [TJ 

53m(G^+*°) = [(l + Axl/fcV]/*G^;^0)"il*53m(G^+^°)(l + Avl/,v^^^^ 
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and (1 + A^fc^^,G^+*°)-i leaves the space spanned by and e^+L invariant. Therefore one also 
obtains S5m(^*G^+*°^fc) > showing E e A^^^k. 

(ii) From (US}, 

G^^ = G^ + - = - A G^v]/,vi/*G^^ , (34) 

and hence ^^^G^^^^ = (1 + Thus Lemma [6](iii) implies 

e*,G^^e, = (1 + A elG'^ek)-' elG^e^ . (35) 

Thus in the hmit e j 0, e^G^+'^Cfc 00 if and only if ^ -A-^. 

(iii) From (!35l) one deduces that the map A 1-^ /^w,Va,A: is integrable in the *-weak topology 
over intervals [a, 6]. Taking imaginary parts of (!35|) . one obtains 



(1 + A3fJe(e*G'^efc))2 + (A 53m(e*G'f,efc))2 

Let X = 3fte(e^G'yefc) and ?/ = ^m{elGyek) ■ Then arctan( ^ ~'~ f ) anti-derivative of the 
function A t-^ $5m(e^(jy^efc). Therefore J^^dX^m{elGy^ek) is bounded by tt and the integral 
over the whole real line exists and is equal to vr. This means that the integral dX^^y^^k 
actually converges to the Lebesgue measure which has no singular part. 

Now let i? be a set of Lebesgue measure zero. Then dX iicj,Vx,k{.B) = 0. As the measures 
are positive this means that for Lebesgue a.e. A G M one has Hciy^^k^B) = 0. □ 

Note that the equation proved in part (iii) above, dE = dX iii:^y^^k{.dE), is well-known from 
the theory of rank one perturbations. 

Theorem 8 Let uo = (cu, V) he fixed such that the matrices ^>\G^'^'^'^ k, as well as 

^fc^v^*°^fe e3;^s^ and are invertihle for Lebesgue almost all E. Set Va = V + A\I'fc\E'^. Then for 
Lebesgue almost all A G M, the measure ficj.v^-k ^■^ absolutely continuous w.r.t. fi^j.Vx-^- 

Proof. By the Radon-Nikodym theorem we can decompose the measure /U^,VA,fc = fx f^d},Vx,i + /^a 
where fx is a function and fix is the part of fJ'w,Vx,k which is singular to ficd,Vx,i- The statement of 
the theorem is that fix = for Lebesgue almost all A. 

In order to show this, we first need to verify a few identities. By multiplying fl34|) with 
from the left and from the right, one obtains 

nc^vx^i = ii + mG'v^kr'nd^v'^, = "^l^^Ji (36) 

1 + XelG^ek 

where the last identity follows from Lemma [6]^iii). From flMl) . one also obtains 

nGvx"^! = nG'v^^i - XnO^'^knGl^I^i . (37) 
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Inserting (136|) in (l37j) gives 



= - A ^ " 1 ' . (38) 

1 + A elG^ek 



Furthermore, it follows from fl35l) that 



l + Ae*G^e, = ^g^. (39) 

Now let A C M be the set of all E where the limit Gy~^^^ exists and all four matrices 
^!Gv^'°^i,^iG^+*°^fc and are invertible. By assumption, the set 

A has full Lebesgue measure and thus by Proposition [TST iii) we have /iii>,VA,A: = fJ'Cj,VxM^ 
Lebesgue a.e. A e M. Thus we can restrict the measures to the set A. We consider the absolutely 
continuous and singular part of jJi'^yy^^k (w.r.t. the Lebesgue measure) separately and begin with 
the singular part. Inserting fl39|) into fl38|) and dividing by e^GyCfc gives 

Let E & A. Then taking z = E + le and the limit e | 0, it follows that 

hm ^^^^ = lim — A . 

where the last term exists and is not zero (except for A = 0) by the invertibility assumptions for 
E G A. Since |e|.Gy^*^efc| — > cxd as e J. for a.e. E w.r.t. the singular part of /i(i,,VA,fc since, 
by Lemma int^iii), the matrix on the l.h.s. is a multiple of 1, one obtains 



lim^^^^^ ^ 

40 elG^+'^ek 



for every A 7^ and a.e. E & A w.r.t. the singular part of fic^j^^^klA. This implies that the 
singular part of flx\A vanishes for every A 7^ and thus the singular part of vanishes also for 
Lebesgue a.e. A G M. 

It remains to consider the absolutely continuous part of fix- Multiplying both sides of ( l38l) 
with |1 + Ae^GyCfcp and taking imaginary parts gives 

|l + Ae^G^e,f53m(vI/*G^^vl/i) = \1 + \elG'^ek\'Qm (^^illG'^^,^ - \Qm (^^{G^^^ 

+ \' [53m(e*G^efc)3f?e(vi/*G^vl>,vl>*G^vl/,) _ g?e(e*G^e,)53m(v|/*G^v|/,vl>*G^vl/,)] . (40) 

For 2; G Ui, the r.h.s. of ( l40l) is a second order polynomial in A which we denote by P{z, A). For 
z = E + ze and E E A, it converges as e | to a limiting polynomial P{E + lO, A). As above 
consider 

Av,fc = |e G M exists and 5>m(e^G^+'°efc) > 
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Claim: For E & An Ay^k, P{E + zO, A) cannot vanish identically as polynomial in A. 

Suppose the contrary. Then by considering the constant and the linear term one deduces 
$5m(^*G^+*°^i) = and S^m = 0. Finally the quadratic term then 

gives S5m(e^G^+*°efc)3fJe (^^tG^+*°^fc^^G^+*°^i) =0. As E G Ay.fc, this now implies that one 

also has 3fte (^\G^^'^^k^lG^^'Hi^ = so that ^*G^+*°^fc^*G^+*^i = 0. This is not the 
case for E & A and hence the claim holds. 

Hence for E G A n Ay.fc, P{E + «0, A) 7^ for Lebesgue a.e. A G M. As the set of {E, A) where 
this happens is clearly measurable, Fubini's theorem implies that for Lebesgue a.e. A one has 
P{E + zO, A) 7^ for Lebesgue a.e. E ^ An Av,fc. Since |1 + Ae^Gy+^^Cfc^ exists and is strictly 
positive for any A G M and E & Af\ Av,fc, it follows from (HOl) that for a.e. A G M, Lebesgue a.e. 
E E An v4v,A:; ^'^(eiGy^*°ei) exists, is finite and strictly positive. Therefore for a.e. A G M, 
the absolutely continuous part of /it:,,VA,i has almost surely a positive density on A fl Ay ^- By 
Proposition [TSVi) the set Ava,A: coincides with Ay^k and, as A has full Lebesgue measure, one 
obtains that A fl Av,^ is an essential support of ^cj,v^,k,a.c- Therefore for a.e. A G M, ficj,Vx,k,ac is 
absolutely continuous w.r.t. ficj,VxA,a.c- This means that also the absolutely continuous part of fl\ 
must vanish for a.e. A G R. □ 

Corollary 3 For fixed Co and Lebesgue a.e. V G Jso*{2L), the matrix valued measure is 
absolutely continuous w.r.t. fi^j,!- Hence for P almost all u = {oJ,V) the measure fi^ is absolutely 
continuous w.r.t. 

Proof. Let uj be fixed. By Proposition [121 the assumptions of Theorem [8] are fulfilled for a.e. 

V G Jso*{2L). Therefore for a.e. V G (]R\I^fc^l/fc)"^, the orthogonal complement of ]R\l'fc\I^fc in 
i7so*(2L), there is some A such that Va = V + A\E'fc\E'^ fulfills the assumptions of Theorem [HI 
Theorem M now states, that for a.e. A G M, the measure fJ'i:j,Vx,k is absolutely continuous w.r.t. 
/^(i,VA,i- For fixed u), the map V h-> {ficj,v,k, f'Cj,v,i) is Borehan as is the Lebesgue decomposition 
for finite measures which maps v) to the singular part of /i w.r.t. v. Hence the set of V where 
IJ'uy^k is absolutely continuous w.r.t. /U^,v,i is measurable. Therefore Fubini's theorem now 
implies that this set has full Lebesgue measure on Jso*{2L). This holds for any k = 2, . . . ,L. 
As a finite intersection of sets of full measure is still a set of full measure we obtain that for a.e. 

V G i7so*(2L) the measure J2k=i f"^,v,k is a.e. w.r.t. ficj,v,i, namely ficj,v is a.e. w.r.t. fidi,v,i- 

The maps u; t— > /loj and u 1— >• fi^^i are Borelian. By the same arguments as above the set of 
uj = {lj, V) where /i^; is absolutely continuous w.r.t. fi^^i is measurable. As the distribution 
of V is absolutely continuous, we obtain that for any fixed u, for pv almost every V, ficj,v is a.e. 
w.r.t. yUtj^v,!- By Fubini's theorem, we obtain that this is true for P almost all u. □ 

Theorem 9 For P almost every uj one has yUa;,i,sing(I^) = 0. Together with Corollary [3] this 
implies that for P almost alluj, one has fJ'uj,smg(M) = 0. 

Proof. Let us define A^ = | exists and Tr(^>m(G^+*°)) > 0} as well as A^^k = 

{E I exists and S5m(e^G^+*°efc) > 0}. By Lemma [SJ^iii), one has A^ = [j^^^ A^^k- Clearly 
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is an essential support of the a.c. part of fi^j and A^^k is an essential support of the a.c. part 

of 

By Kotani theory and Corollary [3] for P almost all u the set A^^^ has full Lebesgue measure 
and /i^ is a.c. w.r.t. fi^^i. Take such an = (cj,V). Then as /i^ is a.c. w.r.t. /i^^i the sets 
A^ and A^ i differ only by a set of measure zero and hence M \ A^ i is a set of zero Lebesgue 
measure. Let V be the projection of V orthogonal to 'ifi'^l and py be the distribution of V, 
namely the push forward of Py. Now set Va = V + A\E'i^^ and let Va be a pre-image of Va 
under the Cayley transformation. Then by Proposition [T3] one has for Lebesgue a.e. A G M, 
/x<^^,i(R \ A^^^i) = iJ.uj^,iO^ \ A^^i) = 0, where ujx = ('^,Va). As yUa;;„i,sing(^a.;„i)^ = by the 
definition of A^^^i, this implies /ioj^, i, sing (I^) = 0. Now by Fubini's theorem for P a.e. lj the 
situation described above happens for pv a.e. V. Then for py a.e. V we have /i(i,,VA,i,sing(lR) = 
for Lebesgue a.e. A. Note that Py is absolutely continuous and for fixed u the set of V where 
Ai(i>,v,i,sing(IR) = is measurable, because the map V ^ /^tj,v,i is Borelian as well as the Lebesgue 
decomposition. Fubini's theorem thus implies that for Lebesgue almost every V in the strip 
supp(py) +]R^i\E'^ one has /itj,v,i,sing(I^) = 0. As the distribution of V is supported in this strip, 
this also holds for pv a.e. V. 

As mentioned, this situation happens to be true for P a.e. Co. By the same arguments as 
above the set of uj where /iaj,i,sing(I^) = is measurable. Fubini's theorem now implies that 
/^w,i,sing(K) = for P a.e. u. Since for P a.e. u one also has that is a.c. w.r.t. /i^^^i, we finally 
obtain that fJ^uj,smg{^) = for P a.e. u. □ 

Proof of second claim of Theorem [T](ii). This is Theorem [91 □ 
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